Physical mechanism for a kinetic energy driven zero-bias anomaly in the 

Anderson-Hubbard model 
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The combined effects of strong disorder, strong correlations and hopping in the Anderson-Hubbard 
model have been shown to produce a zero bias anomaly which has an energy scale proportional to the 
hopping and minimal dependence on interaction strength, disorder strength and doping. Disorder- 
induced suppression of the density of states for a purely local interaction is inconsistent with both 
the Efros-Shklovskii Coulomb gap and the Altshuler-Aronov anomaly, and moreover the energy scale 
of this anomaly is inconsistent with the standard energy scales of both weak and strong coupling 
pictures. We demonstrate that a density of states anomaly with similar features arises in an ensemble 
of two-site systems, and we argue that the energy scale t emerges in strongly correlated systems 
with disorder due to the mixing of lower and upper Hubbard orbitals on neighboring sites. 



Introduction. Since the discovery of high tempera- 
ture superconductivity, there has been widespread inter- 
est in the physics of doped Mott insulators. These ma- 
terials, mainly transition metal oxides, have electronic 
properties which are dominated by electron-electron cor- 
relations and which can be tuned by chemical doping. 
Recently, there have been indications that the disor- 
der introduced by doping can have important effects on 
magnetioi, superconducting^, and transport properties^, 
especially near the Mott metal-insulator transition. In 
particular, experiments have found disorder-induced sup- 
pression of the density of states (DOS) at the Fermi 
energy in several transition metal oxides^"—. Inter- 
pretation of these results is generally given in terms of 
Altshuler-Aronov (AA) zero bias anomalies (ZBA), which 
are known to occur in disordered metals^, as well as the 
Efros-Shklovskii (ES) Coulomb gap, which arises in the 
atomic limit. ^ However, while low disorder samples often 
give results consistent with AA, the pattern in strongly 
disordered samples is a failure to fit either of these pic- 
tures. This highlights the importance of developing a new 
picture which captures the physics of strong correlations 
and strong disorder. 

In this work we address the physics of the ZBA in a 
widely used model for strongly correlated systems with 
disorder, the Anderson- Hubbard model (AHM). The 
Hamiltonian is 



H 



Ccj<T + V (e,n, -I- Un,^h,i) , (1) 



where tij = —t for nearest-neighbor sites i and j, and is 
zero otherwise, and Ciu and the annihilation and 

number operators for lattice site i and spin a, and is 
the energy of the orbital at site i. Disorder is introduced 
by choosing from a uniform distribution G [— -^7 ■§■]■ 
The AHM has been extensively studied by numerous 
methods, including Hartree-Fock^°'^^, dynamical mean 
field iheovy^—, and exact method o ^ ' - ^ . However, the 
existence of a suppression of the DOS at the Fermi level 
(a ZBA) in low dimensions was only recently estab- 
lished conclusively via exact diagonalization and quan- 
tum Monte Carlo calculationaiii^. Ref. [l3 found a V- 



shaped suppression of the DOS ai ep- A narrower soft 
gap, of width ^ O.lt, was also found by exact diago- 
nalization of the one-dimensional AHMjiS and was at- 
tributed to long-range correlations However the V- 
shaped anomaly in Ref. Ill with an energy scale ~ t re- 
mains unexplained. This anomaly is distinct from the 
Mott gap, as demonstrated at half filling by its stability 
at increasing disorder strengths after the Mott gap has 
filled in, and by its persistence away from half filling. 

These properties raise a number of questions. First, 
the suppression of the DOS at the Fermi level by a purely 
local interaction is inconsistent with both the ES and 
AA pictures. The ES Coulomb gap comes specifically 
from the long range 1/r interaction, whereas the purely 
local interaction in the AHM produces a flat DOS in the 
atomic limit. The AA anomaly, at the mean- field level, 
comes from exchange terms which vanish for the AHM. 
Second, the energy scale of the AHM anomaly is of order 
t, in contrast both to the Hartree-Fock energy scale U and 
to the strong coupling scale J ixt^ /U . An alternative to 
the existing frameworks is required to gain insight into 
these behaviors. 

In this paper we present a simple physical mechanism 
which results in a ZBA with energy scale t in the AHM: 
the mixing of the lower Hubbard orbital on one site with 
the upper Hubbard orbital on a neighboring site, an effect 
unique to strongly correlated systems. We demonstrate 
that the DOS of an ensemble of two-site systems displays 
a ZBA which shares key features with that found in the 
lattice. By focusing on the subset of the ensemble which 
contributes most strongly to the anomaly for A > [/ near 
half filling, we isolate the underlying physics. The mech- 
anism in this parameter range is confirmed by an analytic 
calculation of the DOS to leading order int/U . We com- 
ment on the behavior outside this parameter range and 
on the implications for experiments. 

Ensemble of molecules In general, disorder reduces the 
importance of kinetic energy relative to interactions, and 
in the limit of strong disorder hopping causes only mini- 
mal changes to the atomic spectrum. Therefore, an obvi- 
ous starting point for a strongly disordered system is the 
atomic limit. Given that the atomic limit in the case of a 
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FIG. 1: (color online) The DOS of the AHM for an ensemble 
of 500,000 two-site systems, t = 1, (7 = 12, A = 20, and 
fj, = U /2 unless otherwise specified, (a)-(d) show the variation 
of the ZBA with (a) t, (b) U, (c) A, and (d) /.i. Displaying the 
full frequency range, (e) and (f) show that at half filling the 
ZBA goes away and is replaced by the Mott gap for [/ > A, 
and below half filling the ZBA gives way to a step in the DOS. 
This step edge is also seen at the lowest doping in (d). 



purely local interaction has no ZBA, a natural next step 
is to consider an ensemble of two-site systems, referred 
to here as molecules. In considering molecules containing 
just two sites, we can expect to capture only the highest 
energy scale of the ZBA but not the low energy behavior 
coming from longer range correlations. We therefore fo- 
cus our attention on the width rather than the depth of 
the resulting anomaly. 

The molecular Hamiltonian is Eqn. (1) with i, j = 1, 2 
only and site energies chosen from a uniform distribution. 
The ensemble-averaged DOS is smooth and, at half fill- 
ing, preserves electron- hole symmetry. Fig. 1 shows the 
resulting DOS for a range of parameters, with panels (a)- 
(d) focusing near the Fermi level. Panel (a) shows a ZBA 
opening with the addition of hopping. Panels (b) and (c) 
show significant variation in the depth of the anomaly 
with interaction strength U and disorder strength A, but 
essentially no variation of the width. In panel (d) the 
anomaly is the same at half-filling (/i — U /2) and below 
(/i = C//4). While the shape of the anomaly is different 
from that seen in the lattice?^ the width of the anomaly 
shares many features with that of the lattice: Most no- 
tably, the width is not only linear in t but roughly equal 
to t. The ZBA is also independent of U , A, and ^. The 
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FIG. 2: (color online) (a) The phase space of a molecule with 
regions labeled according to the number of particles in the 
atomic ground state and according to the alignment of the 
atomic orbitals as shown in (b). 



advantage here is that the system is sufficiently simple 
that we can explore the anomaly's origin. 

Origin of the anomaly The addition of hopping changes 
the atomic DOS of each molecule, and an exhaustive de- 
scription of all these changes is not sufficiently simple to 
be useful. Here we focus on those molecules which make 
the dominant contribution to the anomaly for A > [/ 
near half filling. We begin by describing two useful cri- 
teria for grouping the molecules, and then explore how 
an anomaly develops from one of the resulting groups 
by a mechanism which is unique to strongly correlated 
systems. 

Recalling that the poles in the Green's function oc- 
cur at the transition energies between many-body states 
and that the allowed transitions depend on the number 
of particles in the ground state, a first useful grouping 
of the molecules is in terms of the number of particles 
in the atomic ground state. Fig. [2ja) shows all possible 
values of the site potentials ei and 62 divided into re- 
gions according to the number of particles in the atomic 
ground state. The changes caused by t will be greatest 
when two atomic orbitals are close in energy, and the 
grouping of molecules may be further refined by distin- 
guishing two ways in which this alignment may occur. 
In the atomic limit, interactions split the resonance at 
each site into a lower Hubbard orbital (LHO) with en- 
ergy Ci and an upper Hubbard orbital (UHO) with energy 
ei + U . In a molecule, the LHO of one site may be more 
closely aligned with the LHO of the other site (LHO- 
LHO in Fig.jjfb)) or it may be more closely aligned with 
the UHO of the other site (LHO-UHO in Fig. [^b)).^ 
When A > [/ near half filling, the largest contribution to 
the anomaly comes from molecules with 2-particle atomic 
ground states and LHO-UHO alignment. Fig.[5][a) marks 
these molecules in red, and labels those in which site i 
is doubly occupied in the atomic limit \di) and those in 
which both sites are singly occupied \s). 

Fig. [3ja) shows the contribution to the DOS of 
all molecules in the ensemble except the \d) and \s) 
molecules. While there is a small ZBA, it is dwarfed by 
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FIG. 3: (color online) Contribution to the DOS (a) from 
all molecules except \d) and \s) and (b) from \d) and |s) 
molecules. Inset shows position of maximum in DOS vs. t. (c) 
Analytic expression (Eqn. ([5])) for |d) and |s) contribution. 



that of the \d) and \s) molecules, shown in Fig. |31Jb).The 
linear dependence on hopping of the anomaly width is 
shown in the inset. 

The physical mechanism behind the hard gap for this 
group of molecules is very simple: When the LHO of one 
site is close in energy to the UHO on the neighboring 
site, a single electron placed in the system will largely 
sit on the lower energy site and the 1-particle ground 
state is essentially unchanged by t. However, if a second 
electron is introduced it can spread itself over the two 
sites, lowering the 2-particle ground state energy by an 
amount linear in t. LHO-UHO mixing without mixing of 
all four atomic orbitals is only possible when U/t ^ 1, 
and in this sense this mechanism is specific to strongly 
correlated systems. 

An analytic expression for the DOS contribution of 
these molecules can be developed to leading order in t/U. 
The 16x16 two-site Hamiltonian matrix breaks into di- 
agonal blocks with fixed particle number. Moreover, the 
2-particle block further divides into a triplet block and a 
singlet block. Because the states in the triplet block are 
not coupled by t while those in the singlet block are, the 
2-particle ground state is always a singlet and may be 



ul)' 4t4il)' and 
4,cL)|). When 62 ei + [/ and U » t, 



the contribution of C2^C2||) may be dropped to leading 
order in t/U. The ground state energy is therefore given 
by the lowest eigenvalue of the matrix 



ei + £2 -V2t 
-V^t 2ei + U 



(2) 



When £2 > fJ- > ei + U {\di) in Fig. ^a)), the atomic 
ground state is c|-|.cj^|) and with hopping the ground 

state energy is Eg = 2^~U + 2x — y — \/y^ + ^f^, where 
X = {ei+€2 + U -2n)/2 and y = (e2 - ei - U)/2. When 



y = 0, the ground state energy is lowered by y/2t. 

Equally important to the DOS, the energies of the 1- 
and 3-particle states are not shifted to leading order in 
t/U. At zeroth order in t/U, the DOS near ep of a, single 
molecule in the \di) configuration is 



pH = lF{y)\Siuj~x-Vy^T2^) 



+6{uj -x + Vy^+2fi) )■ (3) 



where ^ V + 3.- + 2,7gTg 

2y2 + 4i2 + 2j/Vy^T2t2 

The contribution of all |c?i) molecules is determined by 
integrating over < ci < fJ. — U and < £2 < y. 
F{y) depends weakly on y and we set it to one, its lim- 
iting value for large y. The two sites in a molecule are 
equivalent, so the contributions of jdi) and \d2) molecules 
are the same. An analogous calculation may be done for 
the |s) molecules. Combining these and normalizing by 
dividing by A^, we obtain the following disorder- averaged 
DOS due to the Id) and Is) molecules. 



where f{uj,A) 



f{uj,Ada)+ f{uj,A 
{\u\+Af-2e 



(5) 



\UJ\ 



2t^ 



and where a = sign uj. When U — A/2 < ji < C//2, 
Ad+ = (A/2) - J7 + ^ and Ad- = (A/2) while As+ = 
As- = U/2. And when {U - /\)/2 < < U ~ A/2, 
there are no molecules of type |d), and Ag^ = U/2 and 
A,_ = (A/2) + /x-(t//2). 

The key point is that p(a;) in Eqn. ([5]) has a peak at 
^/2t. Fig. [3Jc) shows this analytic expression as a func- 
tion of t at half filling for A = 20 and C/ = 12. The linear 
dependence of the peak position on t is clear. Moreover, 
the fact that the position of the peak in the function 
/(w, A) is independent of A means that the width of the 
anomaly is independent of all remaining parameters — C/, 
A and /i — as seen in Fig.[T]and in the lattice calculations. 
At /I = [/ — A/2, two changes occur: First, the atomic 
DOS is no longer flat ep but instead has a positive 
step, and, second, there are no more \d) molecules. These 
effects are reflected in the fi = U/8 curve in Fig. [TJd). 

This calculation retains only the kinetic energy savings 
of the 2-particle ground state, and the similarity of the 
results to those in Fig. [HJb) demonstrates that this is the 
physics underlying the anomaly in the molecular ensem- 
ble. Moreover, the parallels between the behavior of the 
molecular ensemble and that of the lattice suggest that 
this may in fact be the physics underlying the energy 
scale t of the lattice anomaly as well. 

Discussion We have focused here on the \d) and |s) 
molecules which make the dominant contribution to the 
anomaly for A > [/ near half filling. One might ex- 
pect linear t dependence from all molecules by the fol- 
lowing logic: For molecules with LHO-LHO alignment. 
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the 1- and 3-particle atomic states are shifted in energy 
by an amount Unear in t, while the 0-, 2- and 4-particle 
atomic states are unchanged to zeroth order in t/U. Only 
transitions between states which differ in particle number 
by one enter the single-particle DOS, so a linear depen- 
dence on t would seem to arise in all cases. In fact, when 
the contributions of all LHO-LHO aligned molecules are 
summed, there is a cancelation of all terms linear in t. 
As a result, other than the \s) and \d) molecules the only 
molecules making a contribution to the ZBA with lin- 
ear t dependence are those with 1- and 3-particle atomic 
ground states and LHO-UHO alignment. In this case 
the addition of hopping can change not just the transi- 
tion energies but also the number of allowed transitions, 
by causing a change in the number of particles in the 
ground state. A more detailed discussion of the two-site 
ensemble will be provided elsewhere.— 

The range of disorder and interaction strength to which 
this picture is applicable depends on coordination num- 
ber, with higher disorder required in systems with higher 
coordination number. For systems in which the disorder 
strength is much greater than the bandwidth and the in- 
teraction strength is much larger than the hopping, the 
framework presented here suggests the following: First, 
in a strong magnetic field the alignment of electron spins 
would preclude the sharing of the second electron, reduc- 
ing the depth of the anomaly but not its width with in- 
creasing field. Second, a number of experiments^i^ii have 
observed hard gaps. In particular, the energy dependence 
of the anomaly seen at high disorder in Mn doped LaNiO^ 



bears an intriguing resemblance to Fig. [3jb). Particular 
combinations of dopant and parent compound may result 
in more frequent occurance of the LHO-UHO alignment. 
Finally, an anomaly of width t is not expected for [/ ^ A 
where LHO-UHO alignment is not present. In apparent 
contrast, Ref. [13 states that a ZBA persists for U ^ A. 
However, the width of the DOS suppression in the lattice 
when [/ 3> A is no longer t and decreases with increas- 
ing disorder The low energy physics of the anomaly 
and its behavior outside the parameter range to which 
the two-site ensemble is relevant-disorder less than or of 
order the bandwidth-remain important open questions. 
We suggest that at low energies a mechanism analogous 
to the two-site picture may persist: the mixing by hop- 
ping of many-body states, but in this case states which 
extend over multiple sites. 

In summary, we have demonstrated how disorder can 
generate a feature in the DOS of strongly correlated sys- 
tems with an energy scale given by the hopping, rather 
than the usual strong-coupling scale J. The origin of the 
effect is the kinetic energy savings allowed by spreading 
an electron between two neighboring sites. The fact that 
the key mixing occurs between the LHO on one site and 
the UHO on the neighboring site gives this mechanism a 
character unique to strongly correlated systems. 
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